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ABSTRACT 


The  problem  of  allocating  commodity  flows  to  transport  equipment 
units  In  a  mixed  fleet  Is  treated  by  a  set  of  linear  programming  trans¬ 
portation  models.  The  models  allow  for  the  allocation  of  commodity 
flows  In  a  system  which  may  be  regarded  as  one  whole  entity  such  as  a 
firm,  or  as  a  set  of  reaches  differentiated  by  equipment  characteristics 
or  as  an  environment  of  competing  firms,  or  any  combination  of  these 
The  models  were  shown  to  generate  feasible,  optimal  and  integer  solu¬ 
tions  in  their  application  to  a  hypothetical  two  commodity  example. 
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FOREWORD 

This  paper  represents  the  extension  of  a  previous  investigation 
recorded  in  Working  Paper  No.  2:  "Allocation  of  Open  Hopper  Loads  to 
Empty  Covered  Hopper  Returns"  at  the  Pennsylvania  Transportation  and 
Traffic  Safety  Center.  During  his  labor,  the  author  was  supported  by 
the  U.  S.  Army  Corps  of  Engineers,  North  Central  Division,  under 
contract  number  DACW23-72-C-0009. 
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I .  INTRODUCTION 

The  purpose  of  this  study  is  to  investigate  the  problem  of 
allocating  commodity  flows  to  transport  equipment  units  from  a  mixed 
fleet.  Rather  than  encroaching  upon  the  realm  of  scheduling  and  routing 
however,  the  consideration  in  this  paper  is  geared  more  towards  planning 
purposes,  thus  enabling  the  use  of  static  models.  That  is,  the  assump¬ 
tion  inherent  in  the  solution  of  this  allocation  problem  is  that  the 
relevant  supplies  and  demands  for  some  fixed  time  period  are  known  or 
can  be  estimated,  thus  removing  the  dynamic  nature  of  real  time  from 
the  problem.  Bearing  this  assumption  in  mind,  a  set  of  linear  program¬ 
ming  transportation  models  are  presented  with  initial  application  being 
made  on  the  inland  waterway  system. 

The  problem  traces  its  origin  to  the  devel^nmental  efforts  of  a  com¬ 
puter  simulation  model  at  the  Pennsylvania  Transportation  and  Traffic- 
Safety  Center  designed  to  study  the  inland  waterway  transportation 
systems.  Operationally,  this  simulation  model  is  divided  into  two 
parts.  The  first  section  is  a  tow  generation  program  (TOWGEN) ^ ,  which 
produces  a  time  ordered  list  of  tow  departures  throughout  the  system. 

2 

This  list  is  then  processed  by  a  waterway  simulation  program  (WATSIM)*-. 

The  usefulness  of  this  paper  is  not  limited  to  the  simulation 
model  or  the  waterway  systems  in  general  but  in  fact  the  models  presented 
here  are  applicable  to  any  mode  of  transportation.  Since  the  initial 


H?aterway  Simulation  Series.  Vol.  Ill;  TOWGEN:  A  Tow  Generation  Model 
by  Michael  S.  Bronzini. 

2 

Waterway  Simulation  Series,  Vol.  II;  WATSIM:  A  Waterway  Transport  Sim¬ 
ulator  by  John  H.  Gimbel,  III. 
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application  of  these  models  is  made  in  the  waterway  system,  however,  a 
brief  description  of  TOWGEN  is  provided  in  the  next  section  The 
remainder  of  this  paper  is  devoted  to  problem  definition,  a  general 
description  of  the  models,  the  mathematical  formulation,  and  model 
applications. 


i 
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II.  TOWGEN 

The  purpose  of  TOWGEN  is  to  convert  the  input  commodity  0-D  (origin 
to  destination)  tonnage  matrices  into  a  set  of  0-D  movements  of  tows 
having  known  characteristics.  In  doing  so,  TOWGEN  must  exhaust  the  0-D 
tonnage  matrix  for  each  commodity  while  simultaneously  satisfying  the 
balance  principle.  This  balance  principle  is  stated  as  follows: 

"The  numbers  of  towboats  and  barges  of  each  type  which  arrive 

at  and  depart  from  each  port  in  the  system  should  be  equal  in  the 

long  run." 

This  is  simply  a  steady  state  requirement,  and  says  that,  for  each  type 
of  equipment  in  use  on  the  waterway,  input  must  equal  output  at  every 
point. 

Figure  1  presents  a  generalized  logic  diagram  of  TOWGEN.  The  left 
hand  side  of  the  figure  depicts  the  various  computational  operations 
which  occur  within  the  model,  while  the  right  hand  side  shows  the  source 
and  disposition  of  the  data  files  which  TOWGEN  processes.  The  dashed 
lines  connecting  the  two  sides  of  the  figure  represent  information  flow 

The  input  data  required  by  TOWGEN  consist  primarily  of  the 
following: 

(1)  0-D  tonnage  matrix  for  each  commodity; 

(2)  table  of  barge  data,  showing  commodities  carried,  average  load¬ 
ing,  and  dedicated  equipment  percentage  for  each  type  of  barge 
in  use  on  the  waterway; 

(3)  0-D  mileage  table  and  other  system  description  data; 

(4)  frequency  distributions  of  flotilla  size  (number  of  barges)  vs. 
towboat  horsepower  for  each  commodity  or  commodity  group. 
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Figute  1 


TOWGEN  Logical  Flow  Diagram 
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The  0-D  tonnage  matrix,  which  gives  the  number  of  tons  of  each  com¬ 
modity  to  be  shipped  from  each  port  to  every  other  port,  is  perhaps  the 
most  important  input,  and  is  the  prime  determinant  of  the  number  of 
tows  which  will  be  generated  for  processing  in  WATSIM. 

TOWGEN  processes  these  data  in  the  following  manner:  the  Commodity 
0-D  matrix  is  divided  through  by  average  barge  loadings  to  produce  a 
matrix  of  loaded  barge  0-D  movements.  Next,  the  movements  of  empty 
barges  required  to  support  the  loaded  barge  activity  are  determined 
This  is  done  by  (1)  providing  empty  movements  of  dedicated  equipment  and 

(2)  eliminating  any  remaining  imbalances  between  barges  originating  and 
terminating  at  each  port  by  providing  empty  barge  movements  which 
minimize  total  empty  barge-miles  of  travel.  The  total  barge  0-D  matrix 
(loaded  plus  empty  barges)  is  then  divided  through  by  average  flotilla  . 
size  to  determine  the  0-D  tow  movements  required.  This  tow  0-D  matrix 
is  randomly  sampled  without  replacement,  thus  determining  the  origin  and 
destination  ports  of  each  tow  on  the  tow  list  and  the  sequence  in  which 
these  tows  will  be  introduced  into  the  system.  Additional  parameters 
derived  for  each  tow,  via  appropriate  random  sampling  procedures, 
include  the  following: 

(1)  time  of  departure  from  the  origin  port,  determined  by  modeling 
system  wide  tow  departures  as  a  Poisson  process; 

(2)  towboat  horsepower  and  flotilla  size; 

(3)  number  of  loaded  and  empty  barges  of  each  type; 

(4)  net  tonnage. 

In  summary,  it  should  be  noted  that  TOWGEN  itself  is  not  a  simula¬ 
tion  model,  but  rather  is  one  part  of  a  simulation  package.  Starting 
with  the  basic  elements  of  commodity  transportation  demand  and  transport 
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fleet  supply,  TOWGEN  uses  an  array  of  analytical  and  Monte  Carlo  tech¬ 
niques  to  determine  the  waterway  transport  demand  in  terms  of  discreet 
traffic  units  (i.e.,  barge  flotillas). 
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III.  PROBLEM  DEFINITION 


The  problem  considered  in  this  paper  begins  with  the  second  opera¬ 
tion  performed  in  TOWGEN,  namely  the  computation  of  empty  barge  movements 
required  to  support  the  loaded  barge  activity  and  to  satisfy  the 
balance  principle.  TOWGEN  first  considers  dedicated  equipment  movements. 
Any  remaining  imbalances  are  then  treated  as  a  linear  programming  trans¬ 
portation  problem. 

1.  Dedicated  Equipment 

Barges  are  considered  to  be  dedicated  if  they  move  loaded  from  an 
origin  port  to  a  destination  port,  and  are  then  returned  empty  to  the 
origin  port.  For  each  barge  type,  the  percentage  of  loaded  movements 
utilizing  dedicated  equipment  is  specified,  and  TOWGEN  calculates  the 


corresponding  empty  barge  movements  as , 

E.,,  =  L.  ..  D, /100  i 

jib  ljb  b 


i ,  j  •  •  •  *  k 

b  =  1 ,  .  .  .  ,  B 


where 


=  number  of  loaded  type  b  barges  moving  from  port  i  to  port  j  ; 

E  k  =  number  of  empty  type  b  barges  moving  from  port  j  to  port  i; 

=  dedicated  equipment  percentage  for  barge  type  b; 
k  =  total  number  of  ports  in  the  system. 

If  all  D^  were  equal  to  100,  there  would  be  no  barge  input-output 
imbalances  at  any  port  in  the  system.  In  all  other  circumstances,  how¬ 
ever,  there  will  normally  be  further  balance  requirements  to  be  met. 

2.  The  Balance  Problem 

To  determine  what  additional  empty  barge  flows  are  needed,  TOWGEN 
computes  the  total  number  of  barges  of  each  type  which  originate  and 


terminate  at  each  port-  In  the  general  situation,  some  ports  will 
originate  more  barges  than  they  terminate,  and  hence  will  have  a  "demand" 
for  empty  barges.  The  reverse  situation  will  hold  true  for  other  ports, 
which  thus  will  have  a  "supply"  of  empties.  Also,  the  total  supply  and 
demand  will  be  equal  for  e<;'h  barge  type.  Consequently,  a  balance  be¬ 
tween  input  and  output  of  barges  at  each  port  may  be  achieved  by 
solving  the  following  problem  for  each  barge  type:  "What  movements  of 
empty  barges  should  occur  so  as  to  exhaust  all  supplies  and  meet  all 
demands?" 

The  problem  posed  above  is  solved  in  TOWGEN  by  determining  those 
empty  barge  0-D  movements  which  minimize  total  empty  barge-miles  of 
travel  throughout  the  entire  waterway  system,  using  the  linear  program¬ 
ming  transportation  algorithm. 

The  approach  used  in  TOWGEN  neglects,  however,  the  opportunity  to 
reduce  costs  through  barge  intermix.  That  is,  a  more  efficient  utiliza¬ 
tion  of  equipment  is  derived  through  the  allocation  of  loads  to  barges 
which,  in  the  absence  of  available  loads,  would  be  required  to  move 
empty  because  of  imbalances  in  the  movements  of  the  primary  commodity 

carried.  This  idea  was  introduced  in  two  "Working  Papers"  produced  by 
3 

the  project  staff.  Essentially,  this  technique  involves  the  utiliza¬ 
tion  of  partial  crosshaul  operations  to  meet  some  of  the  empty  barge 
requirements  prior  to  balancing  barge  movements  via  the  linear  program¬ 
ming  transportation  algorithm.  That  is,  for  three  ports  A,  B,  and  C, 


3 

These  papers  have  been  included  in  the  Waterway  Systems  Simulation 
series  as  Technical  Memorandum  No.  1. 
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movements  of  loaded  barges  from  A  to  C  result  in  a  supply  of  empty 
4  5 

barges  at  C.  *  These  empties  were  then  allowed  to  move  commodities 
from  C  to  B,  then  travel  empty  from  B  back  to  A.  The  entire  array  of 
such  possibilities  was  formulated  as  a  general  linear  programming 
problem,  with  cost  minimization  as  the  objective. 

This  paper  presents  five  models  designed  to  approach  the  problem 
in  a  systematic  way.  Four  of  the  models  are  actually  modifications  of 
the  first  general  purpose  model,  presented  to  take  into  account  the 
degree  of  optimality  desired  and  to  simulate  various  real  world  situa¬ 
tions.  Essentially,  these  models  have  the  following  properties: 

(1)  they  present  a  general  multiport,  multicommodity,  multi¬ 
equipment  model; 

(2)  they  consider  all  possible  combinations  of  movements  to 
produce  the  optimal  set  while  satisfying  the  supply  and  demand 
characteristics  at  every  port; 

(3)  the  technique  mentioned  above  and  that  used  in  TOWGEN  view 
the  waterway  operation  as  a  set  cf  movements  from  port  i  to 
port  j  to  be  optimized  so  as  to  result  in  minimum  total  cost. 
In  addition  to  that  the  models  presented  in  this  paper  view 
the  waterway  system  from  the  equipment  user's  point  of  view. 

Before  presenting  these  models,  it  may  be  useful  to  restate  the  problem 
in  an  example.  Consider  two  commodities,  for  simplicity's  sake,  one  of 
which  needs  to  be  carried  in  special  equipment,  say  barge  type  1.  The 


4 

B  is  between  A  and  C. 

^Considers  only  one  stop  between  any  two  ports. 


other  commodity  can  be  carried  in  the  lesser  expensive  barge  type  2  and 
if  convenient  in  barge  type  1.  Consideration  of  the  commodity  matrix 
for  only  barge  type  1  results,  using  the  classical  linear  programming 
transportation  algorithm,  in  a  minimum  number  of  loaded  and  empty  barge 
movements.  The  allocation  of  the  second  commodity  loads  to  the  empty 
type  1  barges  when  economical,  leads  to  a  reduction  in  the  number  of 
barge  type  2  movements.  The  question  is:  when  is  it  convenient,  i  e  , 
when  is  it  economical,  to  carry  commodity  type  2  in  type  1  barges? 

The  following  section  presents  the  general  description  and  the 
theoretical  justifications  for  each  model.  As  mentioned  before,  these 
models  are  not  meant  to  be  mutually  exclusive  and  in  fact  it  would  not 
be  unusual  to  combine  these  models  to  suit  a  given  situation. 
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IV.  GENERAL  DESCRIPTION  OF  MODELS 

A.  The  General  Model 

The  most  efficient  operation  results  when  the  system  is  viewed  as 
one  unit.  The  system  consists  of  ports  with  individual  supply  and  demand 
characteristics  resulting  in  a  prior  specification  of  certain  origin- 
destination  commodity  movements.  However,  these  movements  need  not  be 
reproduced  as  originally  specified  -in  be  manipulated  around  so  as 

to  result  in  a  total  cost  minits’ * «;  .  This  total  minimization  of  0-D 
movements  is  valid  when: 

(1)  the  waterway  oper.i:  is'-n  w  being  viewed  from  a  systems  approach 
with  the  primary  ►-  Vfcing  the  reduction  in  congestion  or 
delays  without  regard  for  the  resultant  0-D  movements; 

(2)  or  correspondingly  when  the  waterway  operation  is  under  the 
jurisdiction  of  one  corporate  entity  or  to  use  the  more  common 
term  monopoly. 

Obviously  this  model  is  unacceptable  as  regards  to  its  application 
to  the  real  world  waterway  system.  However,  it  serves  two  useful 
purposes.  First,  it  provides  the  general  framework  from  which  other 
models  with  suitable  modifications  can  be  obtained.  Second,  it  provides 
a  good  basis,  in  fact  the  best  considering  it  is  the  most  efficient,  for 
the  comparison  of  the  performance  of  other  models. 

This  model  will  also  be  known  as  the  total  cooperation  model  in 
view  of  the  fact  that  the  monopolistic  firm  can  cooperate  with  itself 
completely.  The  term  "cooperation"  will  be  seen  to  be  useful  in  de¬ 
scribing  the  other  models  since  it  not  or.ly  refers  to  the  amount  of 
competition  present  in  the  system  (total  cooperation  in  this  context  is 
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equivalent  to  zero  competition) ,  but  also  in  fact  points  out  the  degree 
of  optimality  that  can  be  attained  in  the  solution. 

B.  The  Non-cooperative  Model 

This  model  does  not  refer,  as  the  name  may  imply,  to  perfect  com¬ 
petition  but  rather  to  the  fact  that  the  degree  of  optimality  attained 
is  minimal  compared  to  the  other  models.  The  model  is  non-cooperative 
in  the  sense  that  the  originally  specified  0-D  movements  must  be  pre¬ 
served  in  the  final  solution,  and  the  allocation  of  type  2  commodity 
loads  to  barge  type  1  in  the  illustrative  example  mentioned  before  must 
be  made  under  this  condition.  This  condition  implies  that  we  can  no 
longer  consider  only  the  total  supply  and  demand  characteristics  of 
every  port  but  in  fact  must  also  reproduce  the  original  0-D  movements. 
This  model  is  appropriate  when: 

(1)  the  waterway  system  under  consideration  is  being  used  by  a 
set  of  small  firms,  for  example  private  shippers,  who  are 
primarily  concerned  with  their  own  supply  and  demand  rather 
than  the  supply  and  demand  of  each  port  in  the  system, 

(2)  there  are  extraneous  and  unknown  costs  involved  in  manipulating 
the  originally  specified  0-D  movements  and  when  these  costs 
are,  for  the  purpose  of  this  model,  prohibitive 

This  model  has  limited  use  in  application  to  the  real  world 
phenomenon  but  is  useful,  as  will  be  seen  in  the  next  section,  in 
illustrating  the  type  of  modifications  that  can  be  made  on  the  general 
model. 

The  next  two  models  are  known  as  the  partial  cooperation  models  in¬ 
dicating  that  there  is  a  reduction  in  the  degree  of  optimality  attained 


(compared  to  the  general  model)  as  a  penalty  for  including  some  of  the 
real  world  constraints. 

C.  The  Multi-reach  Model 

Reach  definition:  A  system  may  be  considered  as  a  set  of  reaches 
if  and  only  if  the  degree  of  cooperation  for  intrareach  movements  is 
greater  than  that  for  interreach  movements.  In  particular,  the  intra¬ 
reach  movements  must  be  formulated  as  the  general  model  (cooperation  is 
total)  while  the  interreach  movements  must  be  modeled  as  non-cooperative 

As  an  example,  consider  the  hypothetical  construction  shown  in 
Figure  2.  Suppose  Airline  A  has  the  sole  franchise  for  the  charter 
trips  between  Boston  and  Cape  Cod  while  Airline  B  serves  Cape  Cod  and 
New  York.  Both  airlines  use,  to  a  large  extent,  DC-8s  while  the  Boston- 
New  York  route,  served  by  a  number  of  major  domestic  trunk  carriers 
utilizes  the  B-727s. 

In  the  non-dynamic  context  (i.e.  ,  disregarding  the  time  factor), 
New  York  originates  720  trips  but  terminates  only  570  trips  Thus 
New  York  has  a  demand  for  150  plane  trips  while  similarly  Boston  has 
a  supply  of  130  plane  trips.  The  general  model  cannot  be  applied 
however,  since  these  supplies  and  demands  involve  both  DC-8s  and 
B-727s.  That  is  the  150  plane  trips  originating  from  New  York  must  be 
distinguished  as  to  their  destination  between  Riverdale  and  Cape  Cod 
or  Boston  so  that  the  appropriate  type  of  plane  trips  can  be  considered. 
The  non-cooperative  model  is  also  inappropriate  since  the  two  routes 
Boston  to  Cape  Cod  and  Cape  Cod  to  New  York  can  each  be  modeled  as  the 
general  model.  Thus  it  becomes  convenient  to  model  Boston,  Quincy  and 
Cape  Cod  as  one  reach  and  Cape  Cod,  Riverdale  and  New  York  as  another 
reach.  All  movements  except  the  Boston-New  York  transits  have  now  been 
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Hypothetical  System  Construction 
for  the  Multi- reach  Model 


accounted  for  in  one  of  these  two  reaches.  These  transits,  utilizing 
the  B-727s  would  be  classified  as  interreach  movements  and  must  occur 
as  specified,  requiring  an  extra  backhaul  of  100  B-727  vacant  plane 
trips  from  Boston  to  New  York. 

The  model  is  appropriate  when: 

(1)  the  reaches  are  so  geographically  separated  as  to  make  consid¬ 
eration  of  the  whole  system  as  one  unit  invalid.  Movement 
optimization,  therefore,  is  possible  within  the  reaches  but 
not  between  reaches; 

(2)  the  commodities  being  hauled  within  the  reaches  are  compatible 
in  the  sense  of  making  cooperation  possible  but  any  movement 
between  reaches  involves  a  different  incompatible  commodity; 

(3)  the  equipment  required  for  an  interreach  movement  is  incompat¬ 
ible  (physically  or  financially)  with  the  equipment  used 
within  the  reaches.  This  occurs  when  small  towboats  and 
barges  are  used  within  the  reaches  but  larger  towboats  and 
barges  are  used  for  a  movement  between  the  reaches. 

(4)  Combinations  of  above.  States  (1),  (2)  and  (3)  above  are  not 

« 

mutually  exclusive,  but  this  model  is  usable  for  combinations 
of  these  states  as  long  as  the  contraints  concerning  the  use 
of  different  barges  to  carry  different  commodities  are 
prespecified. 

D.  The  Multi-firm  Model 

This  model  attempts  to  approximate,  perhaps  more  than  any  other 
model,  the  real  world  phenomenon  on  the  waterways  in  that  it  considers 
the  waterway  system  from  the  user's  point  of  view  rather  than  that  of  a 
bystander. 

Consider  a  set  of  firms  operating  on  the  waterway.  It  would  not  be 
unreasonable  to  assume  that  a  firm  seeks  to  minimize  its  cost  of 
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equipment  use  while  satisfying  its  own  supply  and  demand  characteristics. 
Therefore  a  firm  with  certain  originally  specified  load  movements  from 
port  i  to  port  j  would  be  expected  to  arrange  its  equipment  movements 
so  as  to  result  in  minimum  cost. 

The  most  efficient  equipment  movements  result  when  all  the  com¬ 
modities  are  being  carried  by  one  firm  and  that  firm  employs  an  optimal 
level  of  fleet  mix  to  satisfy  the  supply  and  demand  at  each  port.  This 
is  of  course  the  concept  used  in  the  general  model.  Correspondingly, 
a  number  of  small  firms  operating  without  any  interline  cooperation 
would  tend  to  produce  a  rather  inefficient  set  of  equipment  movements 
from  a  systems  point  of  view.  This  is  precisely  the  non-cooperative 
model  so  that  the  two  extreme  points  on  the  "optimality"  scale  are 
available  for  comparison  with  the  other  partial  cooperation  models. 

This  model  therefore,  simulates  the  real  world  case  of  a  few  large 
firms  where  optimization  within  a  firm's  operation  is  possible  with 
side  benefits  in  smaller  delays  and  reduced  congestion  in  addition  to 
the  minimization  of  an  individual  firm's  cost  of  equipment  use 

E.  The  Dedicated  Equipment  Model 

This  model  is  formulated  to  bring  the  dedicated  equipment  concept 
into  the  analysis  of  empty  barge  movements.  Recall  that  the  dedicated 
equipment  concept  establishes  a  minimum  value  for  the  number  of  empty 
barges,  calculated  as  a  specified  percentage  of  loaded  barge  movements. 
That  is , 


Eijk  ■  “V  (pk>/10° 
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where 


=  number  of  empty  type  k  barges  traveling  i  to  j ; 

=  number  of  loaded  type  k  barges  traveling  j  to  i; 
=  dedicated  equipment  percentage  for  barge  type  k; 


0  <  P.  <100. 
—  k  — 


This  model  corresponds  to  standard  practices  within  the  towing  in¬ 
dustry  since,  for  some  commodities,  such  as  liquid  chemicals,  virtually 
all  shipments  utilize  integrated  tows  which  move  loaded  in  one  direction 
and  empty  in  the  other.  For  other  barge  types,  the  dedicated  equipment 
percentage  reflects  the  operations  of  private  and  contract  carriers 
engaged  primarily  in  accommodating  specific  trades,  or  of  common 
carriers  hauling  select  commodities. 

This  model  is  primarily  a  mathematically  simple  modification  of  the 
general  model.  It  is  useful,  however,  in  allowing  the  reader  to  fcnus 
on  the  types  of  modifications  that  can  be  made  on  the  general  model. 

It  is  further  useful  in  providing  a  rather  simple  way  to  account  for  a 
real  world  situation  and  points  out  the  fact  that  many  of  these  modifi¬ 
cations  presented  in  this  paper  can  and  probably  will  be  performed 
simultaneously  on  the  general  model  so  as  to  provide  in  the  final 
analysis  a  more  complex,  and  a  more  realistic  simulation  of  the  water¬ 


ways  . 
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V.  INPUT  DATA 

This  section  describes  the  classes  of  data  required  by  these 
models.  Since  the  section  is  all  inclusive  (i.e  ,  lists  all  data  re¬ 
quirements  for  all  models) ,  the  use  of  a  particular  model  or  a 
particular  combination  of  models  might  actually  require  only  some  cf 
these  data.  The  choice,  however,  is  not  difficult  to  see 
The  models  require  three  basic  classes  of  input  data: 

(1)  system  data; 

(2)  barge  data; 

(3)  commodity  data. 

A.  System  data 

This  class  of  data  specifies  the  size  and  character  of  the  network 
of  ports  on  the  waterway  system  under  consideration  It  also  specifies 
the  cost  functions  used  in  the  mathematical  formulation, 

1.  Ports  and  Reaches 

A  port  is  defined  as  "any  point  on  the  waterway  which  is  an  origin 
and  a  destination  for  the  movement  of  goods."  This  includes  system  end 
points.  A  reach  is  defined  as  "any  set  of  contiguous  ports  "  Reaches 
as  mentioned  before  are  used  solely  to  designate  portions  of  the  system 
among  which  the  distributions  of  equipment  or  commodity  characteristics 
significantly  differ.  The  number  of  ports  and  reaches  in  a  system  are 
limited  only  by  available  computer  hardware. 

2.  Cost  Functions 

These  models  require  as  input  two  cost  tables  (or  matrices)  for 
each  barge  type.  In  these  tables,  the  rows  and  columns  represent  ports, 
and  each  cell  of  the  matrices  contains  a  scalar  value  For  one  table. 
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the  scalar  value  is  C  ,  *  cost  of  moving  one  barge  load  in  barge  type  k 

ijK 

from  port  i  to  port  j;  i  4  j.  For  the  other  table  with  the  same  barge 
type  the  scalar  value  is, 

cj.,  =  cost  of  moving  an  empty  barge  type  k  from  port  i  to  port  j;  i  4-  j. 
ij  k 

Therefore,  if  k  is  the  number  of  barge  types,  a  total  of  2k  cost  tables 
are  required. 

Two  points  need  to  be  considered  here: 

a.  If  the  cost  is  a  linear  function  of  distance,  then  minimizing 
one  in  linear  programming  is  equivalent  to  minimizing  the 
other. 

Hence,  in  this  case,  the  cost  tables  are  easily  computed  by  multiplying 
an  0-D  mileage  table,  which  gives  the  river  miles  between  pair  of  ports 
in  the  system,  by  a  scalar  and  adding  a  scalar.  That  is, 


\  +  Vir  l** 


where 


d_  =  miles  from  port  i  to  port  j  found  in  each  cell  of  the  0-D 
mileage  matrix; 
k  =  number  of  barge  types; 

A^jB^  =  constants  corresponding  to  the  fixed  costs  and  the  marginal 

cost  of  moving  one  barge  load  one  mile  in  barge  type  k; 

C,,.  =  as  defined  above, 
ijk 

A  similar  relationship  holds  for  empty  barges  and  for  all  other  barge 
types  with  the  substitution  of  appropriate  constants  in  the  function 
given  above. 

b.  If  the  concept  of  firms  is  brought  into  the  analysis  then  two 
conditions  can  prevail; 
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i)  the  cost  functions  can  be  assumed  independent  of  the  user 
firm.  In  this  case  no  change  is  made, 
ii)  If  the  above  assumption  does  not  hold  then  some  additional 
data  is  needed.  If  p  firms  are  considered  in  the  system 
then  with  k  barge  types,  2pk  cost  matrices  are  required 
instead  of  the  2k  matrices  mentioned  above.  In  these 
matrices,  the  rows  and  columns  are  again  origin  and 
destination  ports  respectively,  however,  the  cell  values 
are, 

=  cost  of  moving  one  load  in  barge  type  k  from  port  i  to 
port  j  for  firm  p;  i  ^  j. 


Similarly , 


cost  of  moving  an  empty  barge  type  k  from  port  i  to  port 
j  for  firm  p;  i  #  j  . 


B.  Barge  Data 

This  class  of  data  includes  the  following: 

(1)  the  number  of  barge  types  used  to  transport  the  commodities 
in  the  system; 

(2)  specification  of  the  barge  types  that  can  be  used  to  carry 
each  commodity; 

(3)  list  of  dedicated  equipment  percentages;  (this  could  be 
specified  by  barge  type,  commodity,  reach,  firm,  or  in  general 
for  any  0-D  movement) . 


C.  Commodity  Data 

The  commodity  data  provide  information  concerning  the  loads  to  be 
moved  within  the  system.  The  commodity  data  can  be  given  in  barge  loads 


where  one  barge  load  of  commodity  k  refers  to  the  amount  of  commodity  k 
carried  in  one  barge,  or  if  the  commodity  flows  are  given  in  tons,  then 
the  average  loading  factor  (in  tons  per  barge)  must  also  be  specified 
so  that  the  tonnage  data  can  be  converted  into  barge  loads. 

The  commodity  data  can  include: 

1.  Total  Commodity  0-D  Matrix 

The  cells  of  this  matrix  contain  the  total  barge  loads  to  be 
shipped  from  each  port  to  every  other  port.  Origin  ports  are  specified 
by  the  row  number,  destination  ports  by  the  column  number  (i.e.,  the 
number  in  row  i,  column  j  specifies  flow  from  port  i  to  port  j).  If 
there  are  m  commodities  specified  in  the  system,  then  m  0-D  commodity 
matrices  need  to  be  provided.  Figure  3  illustrates  a  commodity  0-D 
matrix  for  a  system  with  n  ports  and  m  commodities  where, 

=  barge  loads  of  commodity  m  to  be  shipped  from  port  i  to 
port  j ;  i  4  j . 

2.  Commodity  0-D  Matrices  for  Firms 

The  introduction  of  the  firm  concept  into  the  system  necessitates 
some  further  data.  If  p  firms  are  to  be  considered  in  the  system,  then 
at  least  p-1  commodity  0-D  matrices  must  be  specified  for  each  commodity 
in  addition  to  the  total  commodity  0-D  matrix  mentioned  above.  Origin 
ports  are  again  specified  by  the  row  number,  destination  ports  by  the 
column  number.  The  cell  values  are: 

■*ijmp  =  kar8e  loads  of  commodity  m  shipped  from  port  i  to  port  j 
by  firm  p;  i  j . 


Notice  that  if  p-1  0-D  matrices  are  specified  for  each  commodity  then 
the  0-D  matrix  for  the  pC^  firm  for  each  commodity  can  be  calculated 
from  the  relationship: 


k 
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Xijm  =  J^ijmr  i  ^  j,  i,j  =  1,  2 . n  ports 

where  X,,  .  is  the  cell  value  from  the  total  commodity  0-D  matrix.  Ulti- 
ijm 

mately,  a  system  with  p  firms  and  m  commodities  would  require  pm 
commodity  0-D  matrices. 
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VI.  MATHEMATICAL  FORMULATION 

The  procedure  used  to  solve  the  problem  when  the  appropriate  model 
is  mathematically  formulated  is  the  linear  programming  transportation 
technique.  In  order  to  mathematically  formulate  the  model,  the  input 
data  including  the  model  parameters  must  be  specified.  Therefore, 
consider  the  seven  port,  two  commodity,  two  barge  types  system  given 
below.  This  system  is  convenient  because  solutions  using  the  TOWGEN 
and  the  "Working  Papers"  approaches  are  available  for  comparison. 

Consider  two  0-D  matrices  for  covered  hopper  barge  loads  and  open 
hopper  barge  loads  (illustrative  matrix  in  Table  1).^  These  two  matrices 
result  from  the  need  for  carrying  certain  commodities  in  special  equip¬ 
ment,  such  as  grain  in  covered  hoppers,  to  prevent  damage  due  to  the 
elements.  Other  commodities,  such  as  steel,  sand  and  gravel,  etc., 
need  not  be  protected  and  therefore,  open  hopper  barges  will  suffice; 
however,  these  commodities  can  be  carried  in  covered  hopper  barges  if 
it  is  economical  to  do  so.  The  cost  matrices  for  covered  hopper  barges, 
loaded  and  empty  and  open  hopper  barges ,  loaded  and  empty  are  also 
given . 

The  data  given  above  is  the  minimum  amount  of  data  that  can  be  used 
by  one  of  the  models.  In  fact,  it  is  sufficient  for  the  general  and  the 
non-cooperative  models.  Additional  data  will  be  specified  for  the 
partial  cooperation  and  the  dedicated  equipment  models  to  facilitate 
the  mathematical  formulation. 


6 


All  tables  are  collected  at  the  end. 
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The  principle  variables  are  formulated  as  follows: 

XCLC^  *  loaded  covered  hopper  barges  carrying  covered  hopper 
loads  from  port  i  to  port  j ; 

XLC^j  »  loaded  covered  hopper  barges  carrying  open  barge  hopper 
loads  from  port  i  to  port  j ; 

XEC^j  -  empty  covered  hopper  barges  moving  from  port  i  to  port  j; 

XLO^  -  loaded  open  hopper  barges  carrying  open  hopper  loads  from 
port  i  to  port  j ; 

XEO„  =  empty  open  hopper  barges  going  from  port  i  to  port  j ; 

C _ ^  =  cost  of  moving  a  loaded  covered  hopper  barge  from  port  i 
to  port  j ; 

=  cost  of  moving  an  empty  covered  hopper  barge  from  port  i 
to  port  j ; 

C„2  =  cost  of  moving  a  loaded  open  hopper  barge  from  port  i  to 
port  j; 

C^2  =  Cost  of  moving  an  empty  open  hopper  barge  from  port  i  to 
port  j. 


A.  The  General  Model 

Since  the  general  model  is  only  concerned  with  the  total  supply 
and  demand  at  every  port,  the  formulation  of  the  linear  programming 
format  ignores  the  originally  specified  individual  movements.  The 


formulation  begins  with  the  objective  function. 
7  7 


7  7 


Minimize  Z  l  C  .  (XCLC  +  XLC  )  +  I  I  C,,_  XLO,. 
i-1  j-1  3  1J  1J  i-1  j-1 


7  7  1  7  7  1 

+  E  E  Cijl  ^ii  +  Z  1  Ci12  XEOij  ’ 

i-1  j-1  131  13  i-1  j-1  13  z  13 


i-1  j- 
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The  supply  constraints  from  the  0-D  commodity  matrix  for  covered  hoppers 
are , 

7 

£  XCLC, .  =  a,  where  a,  =  row  1  total; 

j-i  13  1  1 


7 

I  XCLC.  =  a.  where  a.  =  row  2  total; 
j=l  J 


7 

£  XCLC...  =  a.  where  a_  =  row  7  total, 

j-i  73  7  7 

The  demand  constraints  from  the  same  table  are, 

7 

£  XCLC.,  =  b,  where  b,  =  column  1  total; 
i=l  11  1  1 


1 

£  XCLC  „  =  b.  where  b„  =  column  2  total; 
1=1  1 11  1 


•  i 

) 

i 

7 

£  XCLC  7  =  b.  where  b?  =  column  7  total.  j 

i=l  '  ' 


The  balance  constraints  for  this  commodity  are, 


7 

£  [(XCLC  +  XLC  +  XEC  ) 


(XCLCj  1  +  XLCjL  +  XEC^)]  -  0; 


J 
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7 

Z 

j=l 


[(XCLC2j  +  XLC2j  +  XEC2j)  -  (XCLCj2  +  XLCj2  +  XECj2>] 


0; 


7 

Z 

j=l 


[(XCLC,.  +  XLC_ . 
7j  7j 


+  XEC.,.)  -  (XCLC .  _  +  XLC...  +  XEC..,)] 
7  3  3  7  j7  J  7 


0. 


These  balance  constraints  simply  state  that  crosshaul  movements  cancel 
each  other  (i.e. ,  the  model  looks  only  at  net  input-output  imbalances 
at  each  port). 

The  supply  constraints  from  the  0-D  commodity  matrix  for  open  hoppers 
are , 


7 

Z  (XLO  .  +  XLC. .)  =  aa. 
j=l  ^  1 


where  aa.  =  row  i  total; 

l 


7 

Z  (XLO^  +  XLC2 J  =  aa2; 


7 

Z  (XLO  +  XLC?  )  =  aa  . 
j=l  3  J 


The  demand  constraints  from  this  matrix  are, 


7 

Z  (XLO^  +  XLC^)  =  bb  1  where  bb1  = 
i=l  J'JL 


column  i  total; 


E  (XLO  +  XLC  )  =  bb  ; 
i=l  1 


E  (XLO_7  +  XLCi?)  =  bb?. 
i=l  '  ' 


The  balance  constraints  for  this  commodity  are, 


E  [(XLO^  +  XEO^)  -  (XLO^  +  XEO^))  =  0; 


E  [(XLO^  +  XE02J  -  (XLOj2  +  XEO_.2)]  =  0; 


/ 

E  [(XLO?j  +  XEO?J)  -  (XL0_.  ?  +  XEOj?)]  =  0. 


In  order  to  guarantee  that  in  the  final  solution  the  elements  along  the 
diagonal  line  will  be  zero  (i.e.,  to  make  sure  that  the  model  will  not 
generate  positive  movements  from  a  port  to  itself) ,  a  few  additional 
constraints  are  necessary. 

XCLClj  =  XLO^  =  XECij  =  XLC^  =  XEO^  =  0  for  all  i=j  , 
i»j  =  ■  •  • ,  7. 

In  running  the  model  on  the  computer,  these  variables  could  simply  be 
omitted  from  the  model  thus  eliminating  these  constraints. 
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Non-negativity  constraint, 

all  X, .  >  0 
ij  ~ 

where  X  .  refers  to  all  variables  with  X  as  the  first  alphabet, 
ij 

B.  The  Non-cooperative  Model 

This  model  requires  the  preservation  of  the  originally  specified 

2 

individual  movements.  Hence  there  are  49  constraints  (n  for  a  system 
with  n  ports)  from  the  covered  hopper  barge  load  matrix  of  the  type, 
XCLC^  =  N  for  all  i,j;  i,j  =  1,2,  .  .  ,7. 

The  balance  constraints  again  are. 


7 

£ 

j  =  l 

7 

£ 


[ (XCLC  .  +  XLC  +  XEC.  )  -  (XCLC  ,+  XLC  +  XEC. 
ij  ij  *j  ji  :i  ii 


[(XCLC^  +  XLC2j  +  XEC2 J  -  (XCLCj2  +  XLC . 2  +  XECj2)] 


0; 


0; 


7 

E  [ (XCLC^j  +  XLC^  +  XEC?j)  -  (XCLC^ ?  +  XLCj?  +  XEC^)]  -  0. 


The  balance  constraints  state  that  while  the  individual  loaded  barge 
movements  must  be  preserved  as  originally  specified,  the  allocation  of 
empty  barge  movements  is  undertaken  with  the  balancing  of  supply  and 
demand  at  every  port. 

The  indivudual  movement  constraints  from  the  open  hopper  barge 
load  matrix  result  in  49  equations  of  the  form, 


XL01J  +  XLC 


*  M 
ij  ij 


for  all  i,j;  i,j  ■  1,2,  .  .  ,  ,  7. 
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The  balance  constraints  are, 


7 

E  [(XLO^  +  XE0lj)  -  (XLOjX  +  XEO^)  ]  =  0; 


7 

E  [ (XLO. .  +  XEO  . )  -  (XLO . 0  +  XEO.,)]  =  0; 
23  23  32  32 


7 

E  [(XLO  +  XEO  )  -  (XLO  +  XEO  )]  =  0. 
_i  '3  ’ 3  3 '  3' 


Specification  of  zero  values  for  the  main  diagonal  elements  result  in, 

XCLC .  .  =  XLC. .  =  XEC , .  =  XLO.  .  =  XEO.  .  =  0  for  all  i=j ; 
ij  ij  ij  ij 

i,j  =  1,2,  .  .  7. 

The  non-negativity  constraints  are, 

all  X±  >  0 


where  X^ , 
ij 


refers  to  all  variables  with  X  as  the  first  alphabet. 


The  objective  function  is  of  the  same  form  as  before. 


C.  The  Multi-reach  Model 

This  model  requires  as  additional  data,  the  specification  of 
reaches  and  the  ports  included  in  each  reach.  Let  the  system  be  axbitrari- 
ly  defined  as  consisting  of  two  reaches,  one  including  ports  1,  2,  3 
and  4  and  the  other  including  ports  5,  6  and  7  keeping  in  mind  that  the 
theoretical  justification  for  this  definition  must  actually  come  from  the 
real  world  phenomenon. 
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Since  this  model  allows  the  consideration  of  only  the  supply  and 
demand  characteristics  of  every  port  within  any  reach,  Table  1  can  be 


modified  to  be  Table  2.  The  shaded  areas  in  Table  2  are  shown  in 

Tables  3  and  4  respectively.  This  is  simply  a  mathematical  statement 

of  the  fact  that  the  load  requirement  matrix  for  each  reach  can  be 

replaced  by  one  cell  (or  terminal)  as  far  as  interreach  movements  are 

concerned.  The  supply  constraints  from  the  covered  hopper  barge  load 
matrix-reach  1  are, 


4 

Z  XCLC  .  =  Z  ; 
j=l  3 


4 

Z  XCLC„.  =  Z- ; 
i-i  2:1  2 


4 

Z  XCLC. .  =  Z. 

.  ,  4j  4 
J=1 


The  corresponding  demand  constraints  are. 


E  XCLCU  -  W  ; 
1  =  1 


4 

Z  XCLCi2  =  W2; 


4 

Z  XCLC, .  =  W 
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The  balancing  constraints  are, 


E  [  (XCLC ,  .  +  XLC,  .  +  XEC  .)  -  (XCLC .  .  +  XLC,.  +  XEC..,)  ]  =  0; 
j  =  1  13  3 1  J1  3l 

4 

E  [(XCLC^  +  XLC^  +  XEC^)  -  (XCLCj2  +  XLC^  +  XEC^)]  =  0; 


E  [ (XCLC. .  +  XLC. .  +  XEC,  .)  -  (XCLC.  .  +  XLC..  +  XEC..)]  =0. 
,=1  4j  4j  4j7  j4  j4  j4 


The  supply  constraints  from  the  covered  hopper  barge  load  matrix-reach  2 


E  XCLC  *  Z  ; 
3=5  2  5 


E  XCLC., .  =  Z,. 
j-5  7 


The  corresponding  demand  constraints  are, 


E  XCLC  *  W  ; 
i=5  5  5 


E  XCLC,,  -  W_. 
i-5  x/  ' 
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The  balancing  constraints  are. 


7 

I  [ (XCLC, .  +  XLC, .  +  XECe.)  -  (XCLC . ,  +  XLC..  +  XEC.,)] 
_c  5j  5j  5j  j5  j5  j5 


=  0; 


7 

Z  [  (XCLC-. .  +  XLC., .  +  XEC,.)  -  (XCLC.,  +  XLC  . ,  +  XEC  .  ,)  ]  =0 
j=5  7j  7j  7j  j7  3  7  3  7 

The  individual  movement  constraints  are  of  two  sets: 

XCLC „  =  N_  for  all  i  =  1,2, 3, 4;  j  =  5,6,7 

and  XCLC  „  =  N  for  all  i  =  5,6,7:  j  -  1,2, 3, 4. 

The  balancing  constraints  are, 

4 

Z  [ (XCLCc .  +  XLC, .  +  XEC,.)  -  (XCLC..  +  XLC . ,  +  XEC . ,) ]  =  0; 
5j  5j  5j'  j o  j5  j5 


4 

Z  ((XCLC-.  +  XLC,.  +  XEC,.)  -  (XCLC.,  +  XLC.,  +  XEC.,)]  =  0. 
,  7j  7j  7j'  3  7  3  7  j7 


This  same  procedure  now  has  to  be  applied  to  the  open  hopper  load 
matrices.  The  supply  constraints  for  reach  1  are, 


4 

Z  (XLO^  +  XLC^)  -  ZZ1; 


The  corresponding  demand  constraints  are, 


E  (XLOi5  +  XLCi5)  =  WW5; 
i=5  L"‘ 
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E  (XL0.7  +  XLC_7)  =  WW?1 
i=5 


The  balancing  constraints  are, 


E  [ (XLO.  +  XEOc  )  -  (XLO . c  +  XEO .  _)  ]  =  0; 
j_5  3.i  5  j  35  j5 


7 

E  [(XLO,.  +  XEO,.)  -  (XLO.,  +  XEO,,)}  =  0. 
73  ?3  37  j  7 


The  individual  movement  constraints  are  of  two  sets: 

XLO.  .  +  XLC,.  *•  M..  for  all  i  =  1,2, 3, 4;  j  »  5,6,7 
ij  ij  ij  J  ’ 

and  XLO. .  +  XLC. .  *  M. .  for  all  i  =  5,6,7;  j  =  1,2, 3, 4 
ij  ij  ij 


The  balancing  constraints  are, 


E  [(XLO^  +  XL05j)  -  (XLOj5  +  XEOj5)]  =  0; 


E  [(XL0?j  +  XE0?j )  -  <XL0j7  +  XE0j7)]  =  0. 


J 


Note:  A  major  assumption  in  the  balancing  constraints  for  this  model 

is  that  equipment  used  for  intrareach  movements  cannot  interchangeably 

be  used  for  interreach  movements.  If  this  assumption  is  unnecessary, 

all  the  balancing  constraints  noted  in  this  model  are  to  be  replaced  by 

the  balancing  constraints  found  in  the  general  model. 

Finally,  the  specification  of  zero  values  for  the  main  diagonal 

elements  and  the  non-negativity  constraints  are, 

XCLC,  .  =  XLC..  =  XEC  .  =  XLO. .  =  XEO . .  =  0  for  all  i-i ; 

ij  ij  ij  lj  ij 

i,j  =  1,2,  .  .  7 

and 

all  X. .  >0 
ij  ~ 

where  X^  refers  to  all  variables  with  X  as  the  first  alphabet.  The 
objective  function  is  of  the  same  form  as  before, 

D.  The  Multi-firm  Model 

In  addition  to  the  data  used  for  the  general  model,  the  multi-firm 
model  requires  the  specification  of  the  firms,  the  cost  matrices  for 
the  firms  and  the  commodity  0-D  matrices  for  the  firms.  Consider  a 
system  with  two  major  firms.  The  0-D  commodity  matrices  for  each  firm 
for  each  commodity  (covered  and  open  hopper  barge  loads)  are  given  and 
the  cost  matrices  are  also  assumed  given. 

Since  the  total  0-D  commodity  matrix  is  the  combination  of  the  0-D 
matrices  for  each  firm,  this  merely  involves  two  general  models,  one  for 
each  firm.  Hence,  the  constraints  given  in  the  general  model  need  to 
be  applied  twice.  In  the  first  application,  the  variables  XCLC^ , 

XLC,,,  XEC,.,  XLO,  and  XEO,  are  replaced  by  XCLC,  . ,  XLC  ...  ,  XEC,  . 
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XLO^^  and  XEO^^  respectively,  where  the  third  subscript  indicates 
that  these  variables  are  for  firm  1.  Similarly,  the  second  application 
involves  the  use  of  variables  XCLC^^*  XLC  XE^ij2*  ^°ij2  ancJ 
XEO^^  respectively.  If  the  general  model  has  K  constraints,  a  multi¬ 
firm  model  with  p  firms  would  in  general  have  pK  constraints. 

Let  the  cells  in  the  cost  matrices  be  denoted  as  follows: 

=  cost  for  firm  1  of  moving  one  loaded  covered  hopper  barge 
from  port  i  to  port  j; 

-  cost  for  firm  1  of  moving  one  empty  covered  hopper  from 
port  i  to  port  j ; 

=  cost  f°r  firm  1  cf  moving  one  loaded  open  hopper  from 
port  i  to  port  j ; 

C^2 ^  =  cost  f°r  firm  1  of  moving  one  empty  open  hopper  from 
port  i  to  port  j ; 

=  cost  for  firm  2  of  moving  one  loaded  covered  hopper  from 
port  i  to  port  j ; 

C^^2  =  cost  for  firm  1  of  moving  one  empty  covered  hopper  from 
port  i  to  port  j ; 

Cf^22  =  cost  for  firm  2  of  moving  one  loaded  open  hopper  from 
port  i  to  port  j ; 

Cij22  =  cost  f°r  firm  2  of  moving  one  empty  open  hopper  from 
port  i  to  port  j. 

The  objective  function  then  takes  the  form, 


3fi 


7  7  1 
minimize  Z  Z  [C,  .....  (XCLC,.,  +XLC, +  C,,.,  XEC . . . 

i=1  i=1  ijH  ijl  ijl  iJH  U1 


4  CtJ21  “-Oljl  4  Cij21  XEOiJl  4  Cijl2  <XCLCij2  4  *LC1J2> 


4  Cijl2  ^132  4  Clj22  XLOl32  4  Cy22  “V1 


E,  The  Dedicated  Equipment  Model 

As  mentioned  before,  this  model  is  introduced  only  to  bang  the 
dedicated  equipment  concept  into  the  analysis  Since  this  concept 
merely  establishes  a  lower  limit  on  a  particular  O-D  empty  barge  move¬ 
ment,  this  model  requires  additional  constraints  of  the  type: 

XLC . .  +  XEC..  >  P.  XCLC . . 
ji  ji  -  ljl  ij 

and 


^i  i  Pij2 


XLO 


il 


where 


»  dedicated  equipment  percentage  for  covered  hopper  barge 
(barge  type  1)  movement  from  port  i  to  port  j , 


and 


P  2  =  dedicated  equipment  percentage  for  open  hopper  barge 
(barge  type  2)  movement  from  port  i  to  port  j . 

The  dedicated  equipment  percentage  could  be  specified,  in  addition  to 
barge  type,  by  reach  and  by  firm.  Hence,  the  constraints  listed  above 
are  in  addition  to  all  the  constraints  (supply,  demand,  balance,  etc.) 
obtained  from  the  particular  model  being  used  whether  it  be  general, 
non-cooperative,  multi-reach,  multi-firm  or  a  combination  of  these 


VII.  MODEL  APPLICATIONS 


A.  System  Description 

The  TOWGEN  and  "Working  Papers"  (to  be  denoted  W#2)  techniques  are 
compared  in  this  section  with  the  general  and  non-cooperative  models. 

The  example  used  for  comparison  purposes  is  the  two  commodity  problem 
mentioned  earlier. 

The  commodity  0-D  matrices  used  in  the  problem  are  given  in  Table  5 
and  6.  The  problem  assumes  a  linear  (i.e.,  sequential)  seven  port 
system  and  the  distance  between  ports  matrix  is  given  in  Table  7.  Cost 
is  assumed  to  be  a  linear  function  of  distance  and  the  cost  functions 
used  are  given  below. 

If  C  =  total  cost  of  moving  a  barge  from  port  i  to  port  j; 

FC  ■  fixed  costs; 

VC  *  variable  costs,  given  in  dollars /mile; 

D  *  distance  between  port  i  to  port  j  ; 
then  C  =  FC  +  VC  *  D. 

(1)  Covered  Hoppers  Barges 

Loaded:  C  =  25  +  4D 
Empty:  C  =  5  +  D 

(2)  Open  Hopper  Barges 

Loaded:  C  =  20  +  3.6D 
Empty:  C  **  A  +  ,9D 

The  cost  functions  were  developed  solely  to  facilitate  solving  the 
problem  and  while  they  may  not  truly  represent  the  actual  transport 
costs,  the  relationship  between  the  cost  functions  are  assumed  to 
parallel  the  actual  relationships. 


The  final  barge  movements  generated  by  TOWGEN,  W#2,  general  and  the 
non-cooperative  models  are  shown  in  Tables  8,  9,  10  and  11  respectively. 
Each  of  these  tables  is  completely  balanced,  i.e. ,  the  number  of  covered 
hopper  barges  and  the  number  of  open  hopper  barges  originating  at  each 
point  is  identically  equal  to  the  number  terminating  at  each  point. 

Each  of  the  cells  in  the  loaded  covered  hopper  barge  matrices  include 
barges  carrying  either  covered  hopper  or  open  hopper  loads.  Table  12 
summarizes  the  statistics  for  total  cost,  distance  in  barge-miles  and 
number  of  barge  movements. 

B.  Analysis 

It  is  seen  from  Table  12  that,  in  comparison  with  TOWGEN  total  cost, 
the  general  model  allows  a  54%  reduction  while  the  W#2  and  the  non- 
cooperative  models  constitute  0.6%  and  1%  cost  reductions  respectively . 
The  reason  for  the  low  values  from  the  latter  techniques  is  due  tc  the 
fact  that  the  opportunity  for  improving  empty  barge  allocation  is  small. 
Table  13  indicates  the  availability  of  covered  hopper  empties  for 
carrying  open  hopper  loads.  This  table  is  derived  from  Table  5  by 
taking  the  differences  between  origins  and  terminations  to  yield  the 
supply  and  demand  for  empty  covered  hopper  barges.  Note  that  the  total 
number  of  covered  hopper  empties  in  this  table  is  the  same  number  gener¬ 
ated  by  the  TOWGEN  technique.  This  number  is  less  than  10%  of  the  total 
covered  and  open  barge  load  movements  indicating  that  the  opportunity 
for  improving  empty  barge  allocation  is  small. 

Despite  this  aspect  of  the  problem  certain  trends  in  the  final 
solutions  are  evident.  The  W#2  and  the  general  model  generate  180  more 
covered  hopper  loads  (all  in  comparison  to  TOWGEN  solution)  and  therefore 


180  less  open  hopper  loads,  while  the  non-cooperative  model  derives  only 
150  more  covered  hopper  loads  and  correspondingly  150  less  open  hopper 
loads.  The  latter  however,  allocates  barge  movements  more  efficiently 
than  the  W#2  technique,  resulting  in  a  slightly  lower  cost  ($1,118,640 
to  $1,122,960).  The  non-cooperative  model  generates  only  70  empty 
covered  barges  as  opposed  to  120  barges  by  the  W#2  technique  indicating 
that  more  of  the  210  available  empty  covered  hopper  barges  are  being 
used  to  haul  open  hopper  loads. 

In  summary  the  following  points  are  noted: 

1.  the  general  and  the  non-cooperative  models  generate  movements 
resulting  in  lower  cost  than  either  the  TOWGEN  or  the  W#2 
techniques ; 

2.  the  TOWGEN,  W//2  and  non-cooperative  techniques  preserve 
individual  loaded  movements  and  this  constraint  coupled  with 
the  scarce  availability  of  empty  covered  hopper  barges  limits 
the  opportunity  for  cost  reduction  in  this  problem. 

C.  Further  Applications 

This  paper  has  presented  five  models  each  model  considering  one 
particular  characteristic  of  the  waterway  system.  The  models  have  been 
formulated  however,  so  as  to  allow  the  derivation  of  combinations  of 
these  models.  The  strength  of  this  approach,  therefore,  lies  not  in  any 
individual  model,  but  rather  in  its  ability  to  analyze  a  complex  situa¬ 
tion  by  identifying  each  of  the  principal  characteristics. 

An  application  of  these  models  is  currently  being  contemplated  for 
the  Illinois-Mississippi  Waterway  Subsystem.  The  Illinois  Waterway 
extends  for  a  distance  of  approximately  326  miles  from  the  confluence  of 
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the  Illinois  and  Mississippi  Rivers  about  38  miles  above  St,  Louis, 
Missouri,  to  Chicago,  Illinois.  About  56  miles  of  the  Mississippi 
River,  from  lock  and  dam  (L&D)  25  above  the  Illinois  through  L&D  26  and 
27  below  it  is  also  Included  in  this  study.  A  schematic  diagram  of  this 
waterway  system  is  shown  in  Figure  4.  The  system  is  defined  to  consist 
of  15  ports,  numbered  as  shown  in  the  figure,  a  junction  point,  and  the 
channel  segments  connecting  these  points. 

All  the  necessary  data  will  be  obtained  through  the  North  Central 
Division  (NCD) ,  Army  Corps  of  Engineers.  The  main  limitation  encounter¬ 
ed  in  this  task  is  the  preparation  of  input  cards  for  the  linear 
programming  computer  program.  Towards  this  end,  a  data  post  processor 
program  is  being  written  by  the  author 

The  procedure  undertaken  in  the  application  cf  the  transportation 
models  to  the  Illinois-Mississippi  Waterway  Subsystem  can  be  indicated 
in  tree  diagram  fashion  as  shown  in  Figure  5.  The  tree  diagram  shown 
Is  for  one  commodity  only.  The  commodity  0-D  matrix  is  broken  down  by 
firms  and  then  by  reaches.  Note  that  the  non-cocperat ive  model  is 
applied  in  two  cases:  (1)  the  operation  of  small  firms  among  which 
there  is  little  barge  intermix,  and  (2)  the  interreach  movements  for 
large  firms.  The  general  model  is  applied  to  the  intrareach  movements 
of  the  large  firms  and  the  dedicated  equipment  model  is  applied  at  the 
end  if  dedicated  movement  is  specified.  The  solution  to  the  problem 
provides  efficient  equipment  movements  that  can  be  compared  with 
observed  data  to  indicate  areas  of  improvement. 

Other  possible  applications  include  routing  air  cargo  carriers, 
rail  cars,  trucks,  etc.  The  transportation  models  provide  a  rather 
simple  approach  towards  increasing  the  utilization  of  a  mixed  fleet. 
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2 
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'  L&D  NO.  25 


.  PEORIA  L&D 


7,° 

.LaGRANGE  L&D 


Ig)  JUNCTION 

Xiao  NO.  26 


L&D  NO.  27 


Figure  4.  II linois-Mlssiss ippi  Ten-Lock  Subsystem 


procedure  as  for  firm 


Figure  5.  System  Application  Tree  Diagram  for  Commodity 
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In  addition  to  the  system  characteristics  already  included,  namely, 
multifirm,  multireach,  dedicated  equipment  and  cooperation  (or  lack  of 
it),  other  constraints  on  fleet  operation  such  as  minimum  and/or  maximum 
load  limits,  specified  cargo  mixes,  etc.  can  be  easily  incorporated  into 
the  models. 

D.  Suggested  Modifications 

An  inherent  assumption  in  the  transportation  models  is  that  the 
demand  for  commodities  in  a  given  period  is  independent  of  time.  This 
assumption  enables  a  static  solution  to  the  equipment  allocation  problem. 
The  removal  of  this  assumption  involves  a  considerably  greater  task  in 
that  inventories  of  equipment  and  commodities  would  be  required  at  each 
node  (place  of  loading  and  unloading)  so  that  dynamic  scheduling  and 
routing  can  occur. 
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VIII.  CONCLUSION 

The  problem  of  allocating  commodity  flows  to  transport  equipments 
units  with  a  mixed  fleet  was  addressed  in  this  paper.  A  set  of  linear 
programming  transportation  models  were  derived,  initially  for  applica¬ 
tion  on  the  inland  waterway  system  although  not  restricted  to  it 

The  models  essentially  allow  for  different  degrees  of  optimization 
in  the  solution  depending  on  the  characteristic  of  the  system  being 
analyzed.  The  characteristic  referred  to  is  the  amount  of  interline 
cooperation  among  the  different  entities  operating  in  the  system  and 
this  cooperation  can  further  vary  by  reaches,  which  are  mutually  ex¬ 
clusive  subsets  of  the  system.  Allowance  was  also  made  for  the  one-way 
flow  of  certain  commodities  such  as  grain  and  petroleum  by  including 
the  specification  of  dedicated  equipment. 

The  performance  of  the  models  for  a  two  commodity  problem  was  com¬ 
pared  with  the  solutions  obtained  using  two  previous  techniques  (see 
footnotes  on  pages  1  and  6) .  The  models  were  indeed  found  to  be 
functional  (i.e.,  generate  feasible,  optimal,  and  integer  solutions)  and 
further  applications  to  the  lllinois-Mississippi  Waterway  Subsystem  is 


contemplated. 
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TABLE  4 

REACH  2:  O-D  COMMODITY  SUBMATRIX 
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COVERED  HOPPER  BARGE  LOADS  COMMODITY  O-D  MATRIX 
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TABLE  6 

OPEN  HOPPER  BARGE  LOADS  COMMODITY  O-D  MATRIX 
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TABLE  10 

NON-COOPERATIVE  MODEL  SOLUTION 

(a)  Covered  Hopper  Barge  Movements 

\  Destination 
Origins. 
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TABLE  11 

GENERAL  MODEL  SOLUTION 
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TABLE  11 

GENERAL  MODEL  SOLUTION  (CONTD  ) 
(b)  Open  Hopper  Barge  Movements 
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TABLE  12 

SUMMARY  STATISTICS 

TOVJGEN 

Solution  W#_2 

Non- 

cooperative 

General 

Cost 

Covered,  Loaded 

$  465,200 

544 ,000 

504,150 

264,900 

Covered,  Empty 

10,450 

6,200 

2,750 

1,000 

Open ,  Loaded 

637,160 

566,600 

602,480 

249,800 

Open ,  Empty 

16.680 

6,060 

9,260 

8.460 

Total 

$1,129,490 

1,122.960 

1,118,640 
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Barge  Movements 

Covered,  Loaded 
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Distance  (in  barge  miles) 
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Open,  Empty 
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Total 
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